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Description 

TECHNICAL FIELD 

[0001 ] The present invention relates to the automatics 
and computation techniques and, more particularly, to 
methods of the preliminary image processing for sharp- 
ening and contrast enhancement. 

BACKGROUND ART 

[0002] One of the methods known in the prior art is 
described in the Russian Patent No.201 5561 , published 
on 16,04.91 , Int. Class G06K 9/40. According to this in- 
vention, the image correction is made basing on the 
analysis of the original image at the processing pixel and 
the local average value over some neighborhood of this 
pixel. 

[0003] The method of that patent smoothes the origi- 
nal image, thus producing the smoothed image contain- 
ing the low frequency components presenting the image 
background. Then the smoothed image is subtracted 
from the original one producing the second image con- 
taining high frequency components without back- 
ground, said second image is then emphasized and 
added to the smoothed image. 
[0004] The disadvantage of this method is that it em- 
phasizes not only the useful signal but also the noise 
containing in high frequency image components, thus 
degrading the quality of the enhanced image. 
[0005] The method according to U.S. Patent No. 
5,038,388, published on 06.08.91 , Int. Class G06K 9/40 
smoothes the original image and subtracts the 
smoothed image from the original one thus producing 
the second image containing the high frequency image 
components only. Said second image isthenadaptively 
emphasized so that the scaling factor is then higherthan 
larger is a difference between the processing pixels and 
its neighborhood. The output image is produced by add- 
ingthe adaptively emphasized second image to the orig- 
inal image, thus sharpening the image without the noise 
amplification. 

[0006] The solution described in this patent does not 
provide any noise suppression as this method can only 
emphasize the high frequency image components that 
may contain a noise. 

[0007] Furthermore, the disadvantage of this method 
known in the art is that it fails to improve considerably 
the sharpness of weak edges as such enhancement re- 
quires to emphasize the high frequency image compo- 
nents in the regions where a difference between the 
processing pixel and its neighborhood is comparable to 
the noise level. Therefore the edge enhancement in 
such regions causes the noise emphasis. 
[0008] According to the image processing method 
disclosed in U.S. Patent No. 5,381,490, published on 
10.01 .95, Int. Class G06K 9/40 the largest difference A 
between the processing pixel and its nearest neighbors 



is calculated. Depending on the magnitude of this differ- 
ence, one of the three processing modes is selected: 

the edge enhancement by means of emphasis of 
5 the high frequency image components if A > T-j, 
here T 1 presents the first pre-defined threshold val- 
ue; 

reproduction of the original image, if T 2 < A < T 1 , 
where T 2 stands for the second pre-defined thresh- 
10 old value. 

- the image smoothing to suppress a noise, if A < T 2 . 

[0009] The first disadvantage of this method is that it 
may emphasize the noise selectively if the difference A 

*5 varies around any of the threshold values for some 
neighboring pixels thus transforming small differences 
between neighboring pixels into larger ones by involving 
different processing modes forthese neighboring pixels. 
[0010] Furthermore, this method fails to provide qual- 

20 ity enhancement of images with different noise magni- 
tudes without tuning as the threshold values 1^ and T 2 
are not selected adaptively. 

[001 1] Another approach to the noise suppression in 
images is described in U.S. Patent 5,563,963, published 
25 on 08.10.96, Int. Class G06K 9/40. The method of this 
patent operates by selecting a plurality of groups of 
neighbors for each pixel of the original image, each 
group being a square containing N x N pixels, here N 
varies for example from 2 to 21 . The processing pixel 

30 may be located at any position inside this square group 
of pixels. The least square best fit plane (the planar ap- 
proximation) is usedthen to approximate pixel values in 
each of selected pixel groups and the new value for the 
processing pixel and the goodness of fit are computed 

35 basing on said approximation for each group. 

[0012] The target pixel of an enhanced image is pro- 
duced by the weighted summing of all the newpixel val- 
ues, the weight of a group is then higher than higher is 
the goodness of fit for this group. 

4Q [0013] The disadvantage of this method is that it fails 
to enhance edges as this method provides the noise 
smoothing only. Furthermore, this method requires 
huge computation efforts to build least square approxi- 
mations by hundreds of groups for each of hundreds of 

45 thousands of pixels. 

[0014] A technique for making low contrast parts of a 
digital image more perceptible is disclosed in the Euro- 
pean Patent EP 0 756 247 A1 , published on 29.01.97, 
Int. CI. G06T 5/20. This method decomposes the origi- 

50 nal image into plurality of images in such a way that 
these images belong to different isotropic frequency 
bands. Thereat the image presenting the low frequency 
components of the original image has lower resolution 
(smaller size) than the image presenting the high fre- 

55 quency components. The emphasis processing is car- 
ried out on these images. The degree of emphasis is 
selected to be lower for pixels having values larger than 
a predetermined threshold value than for pixels having 
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small values. The output image is assembled from the 
processed images. 

[0015] The disadvantage of this method is that it em- 
phasizes noise as well as the useful signal. This method 
fails to provide the edge detection and enhancement as 
the isotropic frequency channels are used. 
[0016] The method disclosed in U.S. Patent 
5,739,922, published on 14.04.98, Int. Classes G06K 
9/40, H04N 1/40 operates by splitting of the original 
color image into three isotropic frequency channels: the 
low frequency image components (LF) , the mediu m fre- 
quency components (MF) and the high frequency com- 
ponents (HF). The adaptive emphasis of the HF com- 
ponents and adaptive suppression of the MF compo- 
nents is then carried out, thereat the multipliers for HF 
and MF image components are then higher than higher 
is a correlation between at least two of three basic image 
colors. The enhanced image is obtained by summing LF 
image components with adaptive ty suppressed MF 
components and adaptively emphasized HF image 
components. 

[0017] The image processing method and apparatus 
described in said patent may have the limited applica- 
tion as they are suitable for color images only as the 
correlation between color components is used for car- 
rying out the image processing. 
[0018] Furthermore, the noise suppression according 
to this invention is significantly limited, as the HF image 
components, that also contain a noise, may be empha- 
sized only and the noise suppression in MF image com- 
ponents is limited because no directional splitting of the 
original image is used. 

[0019] The edge detection and enhancement can not 
be obtained by this method as the isotropic frequency 
channels are used. 

[0020] All these disadvantages degrade the quality of 
enhanced images. 

[0021] The most relevant image processing method 
is described in U.S. Patent No. 5,351 ,305, published on 
27,09.94, Int. Class G06K 9/40. According to this patent, 
a plurality of direction ally filtered images is obtained 
from the original image by applying directional filters in 
the frequency domain. The enhanced image is then 
formed by selecting each target pixel either from the di- 
rectionatly filtered image, if the contrast edge is detected 
nearby the processing pixel or from the original image 
otherwise. Thereat the contrast edge is detected nearby 
the processing pixel by generating the standard devia- 
tion image and by producing the eigenvector description 
of this image. The eigenvector length is compared to the 
predetermined threshold value to detect an edge. 
[0022] The target pixel is equal to the corresponding 
pixel of the original image, if the edge was not detected 
nearby. Otherwise, the target pixel is selected from an 
image filtered with the most nearly corresponding direc- 
tion of filtering. 

[0023] While detecting edges, the eigenvector length 
may vary around the threshold value for several adja- 



cent pixels. Thereby the neighboring pixels of the en- 
hanced image are selected from different images (the 
original image and directionally filtered image) thus 
causing the selective noise emphasis. This emphasis 

5 degrades the enhanced image quality. 

[0024] Furthermore, original images may differ in their 
noise levels thus requiring different threshold values. 
This method does not include the adaptive selection of 
the threshold value and therefore may not provide high 

*0 quality processing of images with different noise levels. 
[0025] Provided that the edge is detected nearby, the 
selection of pixels of the enhanced image is made from 
one of plurality of directionally filtered images thus caus- 
ing the complete suppression of all image structu res that 

'5 differ by their direction from the detected edge, notwith- 
standing that those structures can be clearly seen in the 
original image. 

DISCLOSURE OF INVENTION 

20 

[0026] The object of the claimed invention is to pro- 
vide a method for enhancing the image sharpness and 
contrast combined with simultaneous noise suppres- 
sion. 

25 [0027] The said objective is achieved in the image 
processing method, comprising the presentation of the 
original image as a matrix of the discreet picture ele- 
ments (pixels), splitting of the original image into n fre- 
quency channels, each channel presented by an image 

30 matrix of the same size as the original image, the edge 
detection and assembling of the output image from said 
n frequency channels taking the detected edges into ac- 
count, by splitting of the original image into the low fre- 
quency channel and n-1 high frequency channels and 

35 edge detection by computation in each of n-1 selected 
high frequency channels of the correlation between 
processing pixel and its neighboring pixels, followed by 
the comparison of said correlation value with that for the 
corresponding pixels in other high frequency channels 

40 and with the threshold value for this channel. Basing on 
the results of comparison, the weighting coefficients are 
formed for each pixel of each of n-1 high frequency 
channels, and the assembly of the output image is made 
by summing each pixel from the low frequency channel 

45 with all products of the corresponding (by their location 
in the image) pixels of n-1 high frequency channels by 
their weighting coefficients. 

[0028] The said objective is also achieved by selec- 
tion of m of n-1 high frequency channels (2<m<=n-1) in 

so such a way that they differ one from another in the di- 
rection of principal passing only. Therewithal, the 
weighting coefficients for any of pixels of any of m high 
frequency channels is defined basing on the compari- 
son of its correlation value to the threshold value and to 

55 the correlation values of the corresponding (by their lo- 
cation in the image) pixels of other m-1 high frequency 
channels. 

[0029] The said objective is also achieved by compu- 
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tat ion of said correlation as a product of the processing 
pixel value and the weighted sum of its neighbors, said 
weights being anisotropic thereat the direction of this an- 
isotropy corresponds to the direction of principal pass- 
ing of the processing channel. 
[0030] Furthermore, the threshold value for any of n- 
1 high frequency channels is determined by analysis of 
distribution of values, or absolute values, of pixels in this 
channel. The threshold value may be also determined 
by analysis of distribution of values, or absolute values, 
of pixels of the original image. 

[0031] The claimed as the invention image process- 
ing method is characterized by the following features 
that distinguish it from the most relevant method known 
in the prior art: 

1 . while splitting the original image into frequency 
channels, the low frequency channel and n-1 high 
frequency channels are selected; 

2. the value of correlation between the processing 
pixel and its neighboring pixels is used to detect 
edges in the image. More specifically, the value of 
said correlation is compared to the correlation val- 
ues for the corresponding (by their location in the 
image) pixels in other n-1 high frequency channels 
and to the threshold value for this channel The 
weighting coefficients for all pixels of n-1 high fre- 
quency channels are generated as a result of this 
comparison. 

3. the output image is assembled by means of sum- 
ming of each pixel from the low frequency channel 
and all the corresponding (by their location) pixels 
from n-1 high frequency channels multiplied by their 
weighting coefficients. 

[0032] Concerning the first feature it should be noted 
that the extraction of the low frequency channel that is 
not a subject of any processing provides the distortion 
free reproduction of large-scale objects of the original 
image in the enhanced image, as these objects are 
passed through the low frequency channel without any 
processing. Therewithal the energy of noise passing 
through the low frequency channel without suppression 
is low as most of the noise energy is concentrated at 
high frequencies. 

[0033] Furthermore, the extraction of the low frequen- 
cy channel helps reduce the computation efforts needed 
to find correlation values for n-1 high frequency chan- 
nels, as the subtraction of the local average value is re- 
quired to compute correlation. This subtraction is made 
while extracting (subtracting from the original image) the 
low frequency channel. 

[0034] The splitting of the original image into several 
high frequency channels helps improve, compared to 
selection of one ortwo channels, the noise suppression, 
as a noise associated with pixels of those frequency 
channels where edges were not found nearby, is pre- 
vented from contributing to the enhanced image. 



[0035] Concerning the second feature it should be 
noted that edge detection by means of the correlation 
value between the processing pixel and its neighbors 
allows to find weak edges on a noisy background as the 

5 correlation value is close to zero for the typical noise, 
thereatthe correlation values for adjacent pixels forming 
the edge are positive and of similar magnitude. This dif- 
ference in correlation behavior helps achieve the high 
noise suppression combined with edge emphasis thus 

'0 increasing the image sharpness and contrast. 

[0036] The edge detection basing on the correlation 
between a processing pixel and its neighbors provides 
a quantitative measure of the edge intensity for each im- 
age pixel. This allows performing the selective empha- 

'5 sis of edges having certain intensity, for example, the 
weak edges. Furthermore, use of the correlation be- 
tween the pixel and its neighbors makes the method of 
the present invention applicable to various types of im- 
ages, including, for example, color and gray scale im- 

20 ages, as well as three-dimensional images. 

[0037] The third specific feature, namely assembling 
of the enhanced image by weighted su mming of images 
from all frequency channels, helps to remove complete- 
ly the effect of the selective noise emphasis. This effect 

25 is caused by selection of neighboring pixels from differ- 
ent frequency channels. 

[0038] Furthermore, the determination of the thresh- 
old value by means of analysis of statistical distribution 
of pixels provides a high quality processing of images 
30 with materially different noise magnitudes without 
changing of parameters. 

[0039] Furthermore, use of anisotropic frequency 
channels and anisotropic weights makes the image 
processing method highly sensitive to weak edges. 

35 

BRIEF DESCRIPTION OF DRAWINGS 

[0040] The claimed invention is illustrated by the 
drawings of the apparatus that implements the claimed 
to method. 

[0041] Fig.1 is a block diagram of the apparatus. 
[0042] The preferred embodiment of the sub-units of 
said apparatus are shown in more details at Fig. 2 - 5: 

45 Fig.2 is a block diagram of the frequency channel 
splitting unit (splitting unit). 

Fig.3 is a block diagram of the unit for computation 
of correlations (correlation unit). 

50 

Fig. 4 is a block diagram of the unit for forming of the 
weighting coefficients (weighting unit). 

Fig. 5 is a block diagram of the unit to assemble the 
55 output image (assembling unit). 

Fig. 6 shows an example of the pre-defined channel 
selection matrixes. 
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Fig. 7 illustrates the operation of the frequency 
channel splitting unit. 

Fig. 8 is a graph showing an example of the depend- 
ence of the weighting coefficient on the correlation 
value. 

[0043] Referring to Fig.1 , the apparatus contains the 
image source 1 , the output of said image source con- 
nected to the input of the splitting unit 2. The low fre- 
quency output 7 of said unit 2 is connected to the input 
9 of the assembling unit 5, all other outputs of said split- 
ting unit are connected to the corresponding inputs of 
the correlation unit 3, said other outputs are also con- 
nected to the inputs 10 1 - 10 4 of the assembling unit 5. 
The outputs of said correlation unit 3 are connected to 
the corresponding inputs of the weighting unit 4, the out- 
puts of said weighting unit are the inputs 8 1 - 8 4 of the 
assembling unit 5. Thereat the output of said assem- 
bling unit5 is connected to the input of the memory unit 
6, the output of said memory unit is an output of the ap- 
paratus. 

[0044] Fig. 2 shows more in details the preferred em- 
bodiment of the splitting unit 2. It includes the direct Fou- 
rier processor 1 1 , the input of said processor is connect- 
ed to the output of image source 1, the output of said 
direct Fourier processor 11 is connected to the first in- 
puts of matrix multiplicators 12 0 - 12 4 The second inputs 
of said multiplicators are connected to the correspond- 
ing memory units, said memory units hold the pre-de- 
fined channel selection matrixes illustrated at Fig. 6. Any 
of the matrix multiplicators 1 2 0 - 1 2 4 performs the ele- 
meat-by-element multiplication of matrixes supplied to 
its two inputs. The outputs of said matrix multiplicators 
1 2 0 - 1 2 4 are connected to the inputs of the inverse Fou- 
rier processors 1 4 0 - 1 4 4 . The output of the inverse Fou- 
rier processor 1 4 0 is connected to the input 9 of the as- 
sembling unit 5, and outputs of the inverse Fourier proc- 
essors 14 1 - 14 4 are connected to the inputs 10-, - 10 4 
of the assembling unit as well as to the corresponding 
inputs of the correlation unit 3. 
[0045] Fig.3 shows more in details the correlation unit 
3. The memory unit 15 holds the image of the processing 
frequency channel. The input of said memory unit is the 
input of the correlation unit. It is connected also to the 
input of the noise level measuring unit 20, the output of 
said unit is connected to the first input of divider 19. The 
second input of said divider is connected to the output 
of multiplicator 1 8. The first input of said multiplicator is 
connected to the first output of the memory unit 15. Oth- 
er outputs of said memory unit are connected to the 
weighting summator 17, the output of said weighting 
summator is the second input of multiplicator 18. The 
address input of the memory unit 1 5 is connected to the 
address generator 16. The output of divider 19 is the 
output of the correlation unit. 

[0046] The noise level measuring unit 20 may be im- 
plemented according to the U.S. Patent No. 5,657,401 , 



published on 1 2.08.97, Int. Class G06K 9/40. 
[0047] All the memory units are of the random access 
memory type and they are well known in the art. 
[0048] The weighting summator 17 may be imple- 

5 mented as eight scalar multiplicators (the number of 
said multiplicators is equal to the number of adjacent 
pixels). Any of said scalar multiplicators has two equal 
inputs and one output. The outputs of all said multipli- 
cators are connected to the inputs of summator, the out- 

'0 put of said summator is the output of unit 1 7. The first 
inputs of said multiplicators are the inputs of the unit 17 
and the pre-defined weighting coefficients are supplied 
to the second inputs of said multiplicators. 
[0049] Fig. 4 shows the weighting unit 4. The four in- 
to puts of the weighting unit 4 are the inputs of four round- 
ing units 23 t - 23 4 , the outputs of said rounding units 
are connected to the inputs of the address assembling 
unit 24. The output of said unit 24 is connected to inputs 
of four look-up tables 25^ - 25 4 . The look-up table is a 

20 memory unit that stores the values of the weighting co- 
efficient for any set of the four input correlation values. 
The outputs of the look-up tables 25^ - 25 4 are connect- 
ed to the inputs of the memory units 26-j -26 4 , said mem- 
ory units accumulate values of weighting coefficients. 

25 The address inputs of the memory units 26-j - 26 4 are 
connected to the address generator 27 and the outputs 
of said memory units are connected to the inputs of sum- 
mators 28 1 - 28 4 for averaging of weighting coefficients. 
The outputs of said summators are the outputs of the 

30 weighting unit 4. 

[0050] Fig. 5 shows the assembling unit 5. It consists 
of four multiplicators 29 ^ - 29 4 and summator 30. The 
first inputs 8 1 - 8 4 of said multiplicators are connected 
to the outputs of the weighting unit 4 and the second 

55 inputs 1 0 1 - 10 4 of said multiplicators are connected to 
the outputs of the splitting unit 2. The outputs of multi- 
plicators 29 1 - 29 4 are connected to the corresponding 
inputs of the summator 30 and the input 9 of said sum- 
mator is connected to the output 7 of the splitting unit 2. 

40 The output of summator 30 is connected to the input of 
the memory unit 6 that accumulates an enhanced im- 
age. 

BEST MODE FOR CARRYING OUT THE INVENTION 

45 

[0051] The apparatus implements the claimed meth- 
od as it is described hereinafter. Referring to Fig.1 . the 
input image is generated by the image source 1. The 
MRI u nit may be used, for example, as an image source 

50 I.Said MRI unit produces the image of the cross-section 
of an object, this image being a matrix containing dis- 
creet picture elements (pixels). Said image is carried to 
the input of the splitting unit 2. The operation of the split- 
ting unit 2 is described with a reference to Fig.2and Fig. 

55 7. The input image is transformed to the frequency pres- 
entation by the direct Fourier processor 11. Said fre- 
quency presentation contains the complete information 
about the original image and is represented by the ma- 
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trix of the same size as an input image. This matrix is 
passed to identical matrix multiplicators 1 2 0 - 1 2 4 , these 
matrix multiplicators perform the element-by-element 
multiplication of the frequency presentation of the orig- 
inal image by pre-defined channel selection matrixes. 
Said channel selection matrixes are stored in memory 
units 13 0 - 13 4 . Each channel selection matrix contains 
the multipliers for all spatial frequencies of the frequency 
image presentation. Fig. 6 shows the examples of the 
channel selection matrixes. 

[0052] More specifically, as an image is presented by 
a 2D matrix, its frequency presentation is also a 2D ma- 
trix. Fig. 6a shows schematically a frequency presenta- 
tion matrix. The horizontal and vertical spatial frequen- 
cies vary along axes ^ and ky, respectively. 
[0053] The zero spatial frequency corresponds to the 
constant image density. It is located at the crossing point 
(31) of axes k x and k y . 

[0054] Points 32 and 33 represent the largest spatial 
frequency in horizontal direction. The examples of im- 
ages contributing to these points are shown in drawings 
34 and 35. 

[0055] Similarly, the maximal spatial frequency in the 
vertical direction is located at points 36 and 37; the ex- 
ample of image contributing to these points is illustrated 
by drawing 38. 

[0056] The maximal spatial frequencies are located at 
points 39 - 42. The example of an image contributing to 
these maximal spatial frequencies is shown in drawing 
43. 

[0057] Medium spatial frequency in horizontal direc- 
tion is located at point 45. The example of image con- 
tributing to this point is shown in drawing 44. 
[0058] The location of the spatial frequencies in draw- 
ings Fig.6(b-f) corresponds to the scheme depicted in 
Fig.6a. 

[0059] Fig. 6b shows schematically the pre-defined 
selection matrix for the low frequency channel, said ma- 
trix being stored in the memory unit 13 0 . 
[0060] The dark area 46 is filled by the unit values of 
the matrix elements. This area corresponds to spatial 
frequencies that pass through the low frequency chan- 
nel. The white region is filled by the zero values of the 
matrix elements, therefore the frequencies of the white 
region do not pass through the low frequency channel. 
[0061] Fig.6(c-f) show schematically the selection 
matrixes for four high frequency channels, thereat the 
same notations as in Fig.6b are used. 
[0062] It should be noted that the sum of all channel 
selection matrixes Fig.6(b-f) is the matrix with all ele- 
ments equal to 1 . Therefore all the information from the 
original image passes through at least one channel. 
[0063] Referring to Fig.2, each of the matrix multipli- 
cators 12 0 - 12 4 forms on its output the matrix of the 
corresponding frequency channel inthe frequency pres- 
entation. The inverse Fourier processors 14 0 - 14 4 
transform said matrixes to the coordinate presentation. 
[0064] The direct Fourier processor 11 and inverse 



Fourier processors 14 0 - 14 4 may be implemented bas- 
ing on the Fast Fourier Transform algorithm as de- 
scribed, for example, in: Cooley, J.M., Lewis, RA.W. and 
Welch, RD. (1969) The finite Fourier transform. I.E.E. 

5 E. Trans. Audio Electroacoustics AU-17, 2, 77-86. 
[0065] Fig.7 further illustrates operation of the split- 
ting unit. Fig.7a shows the example of an input image, 
Fig.7(b-f) show the images formed on outputs of the in- 
verse Fourier processors 14 0 - 14 4 , respectively, as a 

10 result of processing of the image shown in Fig. 7a. The 
image of a low frequency channel 7b is carried from the 
output of the Fourier processor 14 0 , being the output 7 
of the splitting unit 2, to the input 9 of the assembly unit. 
The images of four high frequency channels are carried 

*5 from the outputs of the Fourier processors 1 4 n - 14 4 , 
being another outputs of the splitting unit 2, to the cor- 
responding inputs of the correlation unit 3 and to the in- 
puts 10 1 - 10 4 of the assembly unit 5. 
[0066] The further processing of these images will be 

20 described on example of one high frequency channel as 
said processing is identical in all high frequency chan- 
nels. 

[0067] Referring to Fig. 3, the memory unit 15 stores 
the partial image of the processing channel. To compute 

25 the unnormalized correlation value, the center pixel val- 
ue 21 and values of its neighboring pixels 22 are se- 
quentially selected from the memory 1 5. Said values of 
neighboring pixels pass to the input of weighting sum- ' 
mator 17, said summator implements the following op- 

30 e ration on pixel values: 

N 

r = 2>i- X * 

40 

where N is a number of pixels in neighborhood 22 of the 
central pixel (preferably N=8), V f are the pre-defined 
weighting coefficients (preferably V|=1 /8) and X ( are the 
pixel values from neighborhood 22. The muhiplicator 1 8 

45 forms a product of the weighted sum of neighboring pix- 
els and the central pixel value. This product is the un- 
normalized correlation value for the processing pixel. It 
is compared to the threshold value by dividing by said 
threshold value (output of the noise level measuring 

50 means 20) in divider 19. The result of said division is 
compared to 1 .0 in the weighting unit 4. The processing 
described herein is repeated for all pixels of the partial 
image of the processing frequency channel. 
[0068] The image of this frequency channel passes 

55 also to the noise level measuring means 20. The noise 
level from the output of means 20 is used as a threshold 
value to normalize correlation values by divider 19. As 
a result, the matrix containing the correlation values for 
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all pixels of the processing frequency channel is formed 
on output of the correlation unit 3, said correlation values 
being normalized by the threshold valu e for the process- 
ing frequency channel. 

[0069] The correlation values formed by the correla- 
tion unit 3 are carried to the weighting unit 4. Referenc- 
ing to Fig.4, said correlation values for four frequency 
channels pass to inputs of the rounding means 23., - 23 4 . 
Said rounding means decrease the data precision to 4 
or 5 bits. 

[0070] The four rounded values from outputs of 
means 23 1 - 23 4 , each containing 4 or 5 bits, are assem- 
bled into one 16 or 20 bit word by the address assem- 
bling unit 24. The formed address is used as an input 
value forfour look up tables 25 1 - 25 4 . Each of said look 
up tables is a memory unit that stores values of the 
weighting coefficients for any combination of four corre- 
lation values in four frequency channels, thereat such 
combination defines the address formed by the unit 24 
in a unique way. 

[0071] Fig. 8 shows the preferred dependence of the 
weighting coefficient W, in any of the frequency chan- 
nels on the correlation value C| in this channel and cor- 
relation values in other three channels, where A stands 
for the threshold value for this frequency channel. The 
weight Wj depends on the correlation value C, and max- 
imal correlation value L in other three frequency chan- 
nels. This dependence is illustrated at Fig. 8 by the 
curves: 

curve A for C, £ 0.7 L, 
curve B for C, = 0.5 L, 
curve C forCj = 0.3 L, 
curve D for Cj = 0.1 Land 
curve E for C, = 0.01 L 

[0072] The memory units 26 1 -26 4 accumulate values 
of the weighting coeff icients generated by look up tables 
25 1 -25 4 . The address generator 27 and summators28 1 
- 28 4 smooth said weighting coefficients in each fre- 
quency channel. Said smoothing is obtained by sum- 
ming in the summator (for example, 28^ the center val- 
ue of the coefficient and its neighboring values being 
sequentially selected from the memory unit (for exam- 
ple, 26.,) by the address generator 27. The smoothed 
values of weighting coefficients formed on outputs of 
summators (for example, 28 t ) pass to the outputs of the 
weighting unit 4. 

[0073] The operation of the assembling unit is de- 
scribed with a reference to Fig. 5. The values of weight- 
ing coefficients for four frequency channels pass from 
the outputs of the unit 4 to the inputs 8 1 - 8 4 of multipli- 
cators 29 ! - 29 4 and the pixel values of the correspond- 
ing frequency channels are carried from outputs of the 
splitting unit 2 to another inputs 10-, - 10 4 of said multi- 
plicators. The products of said pixel values by the cor- 
responding weighting coefficients generated by multipli- 
cators pass to the inputs of the summator 30. Thereto 



the corresponding pixel value of the low frequency chan- 
nel passes to the input 9 of said summator. Thereby the 
summator 30 calculates the value of the target pixel of 
the output image, said pixels are accumulated in the 

5 memory unit 6. 

[0074] The embodiment described herein illustrates 
the method as applied to 2D scalar images. It is under- 
stood however that the claimed method may be applied 
similarly to 3D images. In this case, in the apparatus 

io used to implement the claimed method: the number of 
frequency channel increases, the 3D Fourier proces- 
sors are used instead of 2D ones and the number of 
neighboring pixels (used, for example, to compute a cor- 
relation value) is 26 instead of 8. 

15 [0075] The claimed method may be applied also to 
processing of vector images, particularly the color im- 
ages. Thereat, the 3 components of a vector presenting 
a pixel value may correspond for example to the inten- 
sity of the 3 basic colors for this pixel. In this case, the 

20 scalar operations on pixel values, like Fourier transform 
and summing, are replaced by the corresponding vector 
operations as it is known in the art and the correlation 
is computed as a scalar product of the center pixel value . 
and the weighted vector sum of its neighbors, thereat 

25 the vector summator contains as many scalar summa- 
tors as the number of vector components. 

INDUSTRIAL APPLICABILITY 

30 [0076] The image processing method according to 
the invention has the following advantages. 
[0077] First, use of the correlation between the 
processing pixel and its neighboring pixels helps detect 
weak edges against a noisy background. This feature 

35 provides high noise suppression in conjunction with the 
emphasis of weak edges, thus significantly improving 
the image quality. 

[0078] Furthermore, use of correlation between the 
processing pixel and its neighbors makes the method 

40 applicable to a wide variety of image types, including 
color images, gray scale images and 3D images. 
[0079] Second, the determination of the threshold val- 
ue by analysis of distribution of pixel values provides a 
high quality processing of images with the material dif- 

45 ference in their noise levels without changing of param- 
eters. 

[0080] Third, the distortion-free reproduction of the 
large-scale image structures is achieved due to the sep- 
aration of the low frequency channel. 

50 

Claims 

1. The image processing method, comprising the 
55 presentation of the original image as a matrix of the 
discreet picture elements (pixels), splitting of the 
original image into n frequency channels, each 
channel presented by an image matrix of the same 
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size as the original image, the detection of edges, 
and assembling of the output image from said n fre- 
quency channels taking the detected edges into ac- 5. 
count, CHARACTERIZED BY splitting of the origi- 
nal image into the low frequency channel and n-1 5 
high frequency channels and detection of edges by 
computation in each of n-1 selected high frequency 
channels for each pixel of the correlation between 6. 
processing pixel and its neighboring pixels followed 
by comparison of said correlation value with corre- 
lation values for the corresponding pixels in other 
high frequency channels and with the threshold val- 
ue for this channel; basing on the results of said 
comparison, the weighting coefficients are formed 
for each pixel of each of n-1 high frequency chan- 
nels, and the assembly of the output image is made 
by summing each pixel from the low frequency 
channel with all the corresponding pixels of n-1 high 
frequency channels multiplied by theirweighting co- 
efficients. 



2. The method according to claim 1, CHARACTER- 
IZED BY determination of said weighting coefficient 
for each pixel of each of n-1 high frequency chan- 
nels by comparison of the corresponding correla- 
tion value to the threshold value. 

3. The method according to claim 2, CHARACTER- 
IZED BY the specific dependence of the weighting 
coefficient on the correlation and threshold values: 



15 
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8. 



the weighting coefficient takes a minimal value 
for correlation values that are significantly 
smaller than the threshold value; 
the weighting coefficient smoothly increases 
from its minimal value to its maximal value for 
correlation values that are close to the thresh- 
old value; 

the weighting coefficient takes its maximal val- 
ue for correlation values that are significantly 40 
larger than the threshold value. 

4. The method according to claim 2, CHARACTER- 
IZED BY the specific dependence of the weighting 
coefficient on the correlation and threshold values: 45 

the weighting coefficient takes a minimal value 
for correlation values that are significantly 
smaller than the threshold value; 
the weighting coefficient smoothly increases so 
from its minimal value to its maximal value while 
the correlation value increases to the second 
threshold value, said second threshold value is 
equal to the product of the first threshold value 
by a pre-defined coefficient; 55 
the weighting coefficient smoothly decreases 
from its maximal value to its limit value while 
correlation value is larger than the second 



threshold value. 

The method according to claim 1 , CHARACTER- 
IZED BY m of said n-1 high frequency channels, 
where 2<=m<=n-1 , being different one from another 
in the direction of their principal passing only. 

The method according to claim 5, CHARACTER- 
IZED BY determination of said weighting coefficient 
for each pixel of each of m high frequency channels 
by comparison of the corresponding correlation val- 
ue to the threshold value and to the correlation val- 
ues for corresponding (by their location in the im- 
age) pixels of other m-1 high frequency channels. 

The method according to claim 1, CHARACTER- 
IZED BY the picture element (pixel) being repre- 
sented by a scalar value characterizing, for exam- 
ple, the image intensity at the corresponding pixel. 

The method according to claim 7, CHARACTER- 
IZED BY computation of said correlation value for 
each pixel by multiplication of said pixel value by 
the weighted sum of its neighboring pixels. 

The method according to claims 8 and 5, CHARAC- 
TERIZED BY use of anisotropic weights for compu- 
tation of said weighted sum of the neighboring pix- 
els, thereat the direction of said anisotropy corre- 
sponds to the direction of principal passing for the 
processing frequency channel. 



10. The method according to claim 7, CHARACTER- 
IZED BY determination of said threshold value for 
35 each of said n-1 high frequency channels by anal- 
ysis of distribution of pixel values in an image of the 
processing frequency channel. 



11. The method according to claim 7, CHARACTER- 
IZED BY determination of said threshold value for 
all said frequency channels by analysis of distribu- 
tion of pixel values of the original image. 

12. The method according to claim 1, CHARACTER- 
IZED BY the picture element (pixel) being repre- 
sented by a vector, components of said vector char- 
acterizing, for example, the intensity of the basic 
colors (red, green and blue) at the corresponding 
pixel. 

13. The method according to claim 12, CHARACTER- 
IZED BY computation of said correlation value for 
each pixel as a scalar product of said pixel vector 
by the weighted sum of vectors representing its 
neighboring pixels. 

14. The method according to claims 1 3 and 5, CHAR- 
ACTERIZED BY use of anisotropic weights for 
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computation of said weighted sum of the neighbor- 
ing pixels, thereat the direction of said anisotropy 
corresponds to the direction of principal passing for 
the processing frequency channel. 

5 

15. The method according to claim 12, CHARACTER- 
IZED BY determination of said threshold value for 
each of said n-1 high frequency channels by anal- 
ysis of distribution of absolute values of vectors rep- 
resenting pixels of an image of the processing fre- '0 
quency channel. 

16. The method according to claim 12, CHARACTER- 
IZED BY determination of said threshold values for 

all high frequency channels by analysis of distribu- '5 
tion of absolute values of vectors representing pixel 
values of the original image. 

17. The method according to claim 1, CHARACTER- 
IZED BY smoothing of the correlation values for 20 
several neighboring pixels before computation of 
the weighting coefficients, said smoothing being im- 
plemented at least in one of n-1 high frequency 
channels. 

25 

18. The method according to claim 17, CHARACTER- 
IZED BY non linear transformation of the correlation 
values prior to the smoothing of the correlation val- 
ues, said non linear transformation remains un- 
changed the values that are smaller or close to the 30 
threshold value, and decreases correlation values 
that are significantly larger than the threshold value. 

19. The method according to claim 1, CHARACTER- 
IZED BY smoothing of the weighting coefficients 35 
over several neighboring pixels, said smoothing be- 
ing implemented at least in one of n-1 high frequen- 
cy channels. 

20. The method according to claim 1 CHARACTER- 40 
IZED BY the input image being p dimensional ma- 
trix of picture elements, where p is greater or equal 

to 3. 

21. The method according to claim 1, CHARACTER- 45 
IZED BY use of different threshold values for differ- 
ent parts of the image, said different threshold val- 
ues being used to form the weighting coefficients at 
least in one of n-1 high frequency channels. 

50 

22. The method according to claims 7 and 21 , CHAR- 
ACTERIZED BY determination of said threshold 
values for different parts of the image and different 
high frequency channels by the analysis of distribu- 
tion of pixel values in the corresponding part of the 55 
image of the corresponding frequency channel. 

23. The method according to claims 1 2 and 21 , CHAR- 



ACTERIZED BY determination of said threshold 
values for different parts of the image and different 
frequency channels by the analysis of distribution 
of absolute values of vectors representing pixels in 
the corresponding part of the image of the corre- 
sponding frequency channel, 

Patentansp ruche 

1. Das Verfahren zur Bildverarbeitung, daB auf der 
Darstellung des Ausgangsbildes in Form einer Ma- 
trix der diskreten Bildelemente (Pikseln), auf der 
Aufteilung des Ausgangsbildes in n Frequenzkana- 
le, dabei e nth a It jedes davon eine Matrix derglei- 
chen GroBe wie das Ausgangsbild, auf der Ermitt- 
lung von Bildgrenzen und auf der Formierung des 
Endbildes aus den obengenannten n Frequenzka- 
nale unter Berucksichtigung der ermitterten Gren- 
zen stutzt, gekennzeichnet dadurch, daft bei der 
Aufteilung des Bildes in die Frequenzkanale ein 
Niederfrequenzkanat und n-1 Hochfrequenzkanale 
abgetrennt werden, zur Abtrennung der Grenzen 
wirdinjedemder n-1 Hochfrequenzkanale fur jedes 
Bitdelement die Korrelation zwischen dem verarbei- 
tenden Bildelement und den es umkreisenden Bild- 
elementenberechnet.derberechnete Korrelations- 
wert wird mit der Korrelationswerten fur entspre- 
chende Bildelemente in anderen Frequenzkanalen 
und mit dem Schwellenwert fur dieses Kana! vergli- 
chen, nach den Vergleichsergebnissen wird fur je- 
des Bildelement jedes der n-1 Hochfrequenzkanale 
ein Gewichtskoeffizient gefonnt, und die Formie- 
rung des Endbildes erfolgt durch die Addierung je- 
des der Bildelmente des Niederfrequenzkanals mit 
alien entsprechenden Bildelementen der n-1 Hoch- 
frequenzkanale, murtipliziert mit ihnen entspre- 
chenden Gewichtskoeffizienten. 

2. Das Verfahren nach Anspruch 1 dadurch gekenn- 
zeichnet, daB der Gewichtskoeffizient fur jedes 
Bildelement jedes der n-1 Hochfrequenzkanale 
durch den Vergleich des entsprechenden Korrelati- 
onswertes mit dem Schwellenwert bestimmt wird. 

3. Das Verfahren nach Anspruch 2 dadurch gekenn- 
zeichnet,dafi der Gewichtskoeffizient mit dem Kor- 
relationswert und mit dem Schwellenwert wie folgt 
verbunden ist: 

Gewichtskoeffizient ist minimal, wenn der Kor- 
relationswert wesentlich niedriger als der 
Schwellenwert ist; 

Gewichtskoeffizient steigt monoton von dem 
minimalen bis hin zu dem maximalen Wert an, 
wenn der Korrelationswert mit dem Schwellen- 
wert vergleichbar ist; 

Gewichtskoeffizient ist maximal, wenn der Kor- 
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relationswert wesentlich grofier als der Schwel- 
tenwert ist. 

4. Das Verfahren nach Anspruch 2 dadurch gekenn- 
zelchnet, dafl derGewichtskoeffizientmitdem Kor- 
relationswert und mit dem Schwellenwert wie folgt 
verbunden ist: 

Gewichtskoeffizient ist minimal, wenn der Kor- 
relationswert wesentlich niedriger als der 
Schwellenwert ist; 

Gewichtskoeffizient steigt monoton von dem 
minimalen bis hin zu dem maximalen Wert an, 
wenn der Korreiationswert bis zu zweitem 
Schwellenwert anwachst, der dem Produkt des 
Schwellenwertes mit dem vorbestimmten Ko- 
effizient gleich ist; 

Gewichtskoeffizient nimmt monoton von dem 
maximalen Wert bis hin zu dem Grenzwert ab, 
wenn der Korreiationswert groBer als der zwei- 
te Schwellenwert ist 

5. Das Verfahren nach Anspruch 1 dadurch gekenn- 
zeichnet, dafl m aus n-1 Hochfrequenzkanale, wo 
2 <, m <, n-1, sich nur durch die haupsachliche 
DurchlaBichtung voneinander unterscheiden. 

6. Das Verfahren nach Anspruch 5 dadurch gekenn- 
zeichnet, daB der Gewichtskoeffizient fur jedes 
Bildelement jedes der m Hochfrequenzkanale 
durch den Vergleich des entsprechenden Korrelati- 
onswertes mit dem Schwellenwert und mit den Kor- 
relationswerten fur die nach der Lage entsprechen- 
de Bildelemente in den anderen m-1 Hochfre- 
quenzkanale bestimmt wird. 

7. Das Verfahren nach Anspruch 1 dadurch gekenn- 
zeichnet, daB das Bildelement (Piksel) eine Zahl 
darstellt, die z.B. die Scharfe des entsprechenden 
Bildpunktes kennzeichnet. 

8. Das Verfahren nach Anspruch 7 dadurch gekenn- 
zeichnet, dan der Korreiationswert fur jedes Bild- 
element durch die Multiplikation dieses Bildelemen- 
tes mit der gewogenen Summe der benachbarten 
Bildelemente berechnet wird. 

9. Das Verfahren nach Anspruchen 8 und 5 dadurch 
gekennzeichnet, da(3 bei der Berechnung der ge- 
wogenen Summe der benachbarten Bildelemente 
anisotrope Gewichte benutzt werden, dabei ent- 
spricht die Richtung der Anisotropie der hauptsach- 
lichen Richtung des verarbeitenden Frequenzka- 
nats. 

10. Das Verfahren nach Anspruch 7 dadurch gekenn- 
zeichnet, daB der Schwellenwert fur jedes der n-1 
Hochfrequenzkanale durch die Auswertung der Zu- 



ordnungen von Werte der Bildelemente dieses Fre- 
quenzkanals bestimmt wird. 

1 1 . Das Verfahren nach Anspruch 7 dadurch gekenn- 
5 zeichnet, daB die Schwellenwerte fur alle Fre- 
quenzkanale durch die Auswertung der Zuordnun- 
gen von Werte der Bildelemente des Ausgangsbil- 
des bestimmt wird. 

10 12. Das Verfahren nach Anspruch 1 dadurch gekenn- 
zeichnet, daB das Bildelement ein Vektor darstellt, 
dessen Komponente z.B. die Helligkert der drei 
Hauptfarben (rot, grtin, blau) kennzeichnen. 

'5 13. Das Verfahren nach Anspruch 12 dadurch ge- 
kennzeichnet, daB der Korreiationswert fur jedes 
Bildelement durch die skalare Multiplikation dieses 
Bildelementes mit der gewogenen Summe der be- 
nachbarten Bildelemente berechnet wird. 
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14. Das Verfahren nach Anspruchen 13 und 5 dadurch 
gekennzeichnet, daB bei der Berechnung der ge- 
wogenen Summe der benachbarten Bildelemente 
anisotrope Gewichte benutzt werden, dabei ent- 
spricht die Richtung der Anisotropie der hauptsach- 
lichen Durchla3richtung des verarbeitenden Fre- 
quenzkanals. 

15. Das Verfahren nach Anspruch 12 dadurch ge- 
kennzeichnet, daB der Schwellenwert fur jedes der 
n-1 Hochfrequenzkanale durch die Auswertung der 
Zuordnungen von absoluten VektorengroBen, die 
die Bildelemente dieses Frequenzkanals darstel- 
len, bestimmt wird. 

16. Das Verfahren nach Anspruch 12 dadurch ge- 
kennzeichnet, daB der Schwellenwert fur alle 
Hochfrequenzkanale durch die Auswertung der Zu- 
ordnungen von absoluten VektorengroBen, die die 
Bildelemente dieses Frequenzkanals darstellen, 
bestimmt wird. 

17. Das Verfahren nach Anspruch 1 dadurch gekenn- 
zeichnet, daB vor der Berechnung der Gewichts- 
koeff izienten mindestens in einem der n-1 Hochfre- 
quenzkanale die Ausgleichung der Korrelations- 
werte fur einige benachbarte Bildelemente durch- 
gefuhrt wird. 



so 18. Das Verfahren nach Anspruch 17 dadurch ge- 
kennzeichnet, daB vor der Mittetung der Korrelati- 
onswerte die unlineare Umformung der Korrelatio- 
nen auf solche Weise durchgefuhrt wird, daB die 
Werte, die niedriger als Schwellenwert oder ver- 

55 gleichbar mit dem Schwellenwert sind, unverandert 
bleiben, und die Werte, die wesentlich groBer als 
Schwellenwert sind, vermindert werden. 
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19. DasVerfahrennachAnspruchl dadurch gekenn- 
zeichnet, daB bei der Formierung des Gewichtsko- 
effizientes mindestens in einem der n-1 Hochfre- 
quenzkanale die Ausgleichung dieses Gewichtsko- 
effizientes nach einigen benachbarten Bildelemen- 5 

ten durchgeftihrt wird. 2. 

20. Das Verfahren nach Anspruch 1 dadurch gekenn- 
zeichnet, daft das Ausgangsbild eine Matrix der 
Bildelementen der Dimension p darstelrt, wo p £ 3. w 

21. Das Verfahren nach Anspruch 1 dadurch gekenn- 3. 
zefchnet, dafi fur die Formierung des Gewichtsko- 
effizientes mindestens in einem der n-1 Hochfre- 
quenzkanale die unterschiedliche Schwellenwerte '5 

fur unterschiedliche Bildteile benutzt werden. 

22. Das Verfahren nach Anspriichen 7 und 21 dadurch 
gekennzefchnet t daB Schwellenwerte fur unter- 
schiedliche Bildteile und fur verschiedene Fre- 20 
quenzkanale durch die Auswertung der Zuordnun- 
gen von Werten derentsprechenden Bildteilen des 
entprechenden Frequenzkanals bestimmt werden. 

23. Das Verfahren nach Anspriichen 12 und 21 da- 25 
durch gekennzefchnet, dan die Schwellenwerte 

fur unterschiedliche Bildteile und fur verschiedene 
Frequenzkanale durch die Auswertung der Zuord- 4. 
nungen von absoluten VektorgroBen der entspre- 
chenden Bildteilen des entprechenden Frequenz- 50 
kanals bestimmt werden. 



Revendlcations 

35 

1. Le procede du traitement de I'image fonde sur la 
presentation de I'image premiere com me une ma- 
trice des pixels discrets; le partage de I'image pre- 
miere en n voies de frequence, dont chacune con- 
tient la matrice de m§me dimension que I'image 
premiere; le degagement des frontieres de I'image 
et la formation de I'image de sortie de n voies de 
frequence indiquees ci-dessus en tenant compte 
des frontieres degagees, se dlstlnguant par le fait, 
que pendant le partage de I'image en voies de fre- 45 
quence on degage une voie de basse frequence et 
(n-1) voies de haute frequence, pour le degage- 5. 
ment des frontieres dans chaque (n-1) voies de 
haute frequence pour chaque pixel on calcule lava- 
leurde fa correlation entrele pixel traite et les pixels 50 
qui I'entourent, compare la valeur de la correlation 
obtenue avec les valeurs de la correlation pour les 6. 
pixels correspondants dans les autres voies de fre- 
quence et avec la valeur seuil pour la voie donnee, 
selon les resultats de lacomparaison forme le coef- 55 
ficient de ponderation pour chaque pixel de chaque 
de (n-1 ) voies de haute frequence et forme I'image 
de sortie a I'aide de la so mm at ion de chaque pixel 



de la voie de basse frequence avec tous les pixels 
correspondants de (n-1 ) voies de haute frequence 
mu hip lies aux coefficients de ponderation corres- 
pondants. 

Le procede selon p. 1 se dlstlnguant par le fait, 
qu'on determine le coefficient de ponderation pour 
chaque pixel de chaque de (n-1 ) voies de haute fre- 
quence a I'aide de la comparaison de lavaleur cor- 
respondante de la correlation avec la valeur seuil. 

Le procede selon p. 2 se dlstlnguant par le fait, 
que le coefficient de ponderation est lie a la valeur 
de la correlation et la valeur seuil de la maniere sui- 
vante: 

le coefficient de ponderation est minimum, si la 
valeur de la correlation est essentielement plus 
basse que la valeur seuil; 
le coefficient de ponderation augmente d'une 
maniere monotone de la valeur minimum jus- 
qu'a la valeur maximale, si la valeur de la cor- 
relation est comparable avec la valeur seuil; 
le coefficient de ponderation est maximale, si 
la valeur de la correlation est essentielement 
plus elevee que la valeur seuil 

Le procede selon p. 2 se dlstlnguant par le fait, 
que le coefficient de ponderation est lie a la valeur 
de la correlation et lavaleur seuil de la maniere sui- 
vante: 

le coefficient de ponderation est minimum, si la 
valeurde la correlation est essentielement plus 
basse que la valeur seuil; 
le coefficient de ponderation augmente d'une 
maniere monotone de lavaleur minimum jus- 
qu'a la valeur maximale, si ta valeur de la cor- 
relation augmente jusqu'a la deuxieme valeur 
seuil egale au produit de la valeur seuil et du 
coefficient predetermine; 
le coefficient de ponderation diminue d'une ma- 
niere monotone de la valeur maximale jusqu'a 
la valeur seuil, si la valeur de la correlation est 
plus elevee que ta deuxieme valeur seuil. 

Le procede selon p. 1 se dlstlnguant par le fait, 
que m de (n-1) voies de haute frequence, ou 2 <, m 
^ n-1, ne se distinguent que par la direction de la 
transmission principale. 

Le procede selon p. 5 se dlstlnguant par le fait, 
qu'on determine le coefficient de ponderation pour 
chaque pixel de chaque m voies de haute frequen- 
ce par la comparaison de la valeur correspondante 
de la correlation avec la valeur seuil et les valeurs 
de la correlation des pixels correspondant par les 
positions de m-1 voies de frequence restantes. 
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7. Le precede selon p. 1 se dlstlnguant par te fait, 
que le pixel represents un seul nombre, qui carac- 
terise, par exemple, la brilliance du point correspon- 
dant de f image. 

8. Le procede selon p. 7 se distlnguant par le fait, 
qu'on calcu le la valeur de la cor re I at io n pou r chaque 
pixel par la multiplication de ce pixel par la somme 
ponderee des pixels voisins. 

9. Le procede selon p. 8 et p. 5 se dlstlnguant par le 
fait, qu'on utilise des poids anisotropiques pendant 
la calculation de ta somme ponderee des pixels voi- 
sins, la direction de fanisotropie correspondant a la 
direction de la transmission principale de la voie de 
frequence traitee. 

10. Le procede selon p. 7 se distlnguant par le fait, 
qu'on determine la valeur seuil pour chaque de (n- 
1 ) voies de haute frequence a I'aide de ('analyse de 
la repartition des valeurs des pixels de cette voie 
de frequence. 

11. Le procede selon p. 7 se distlnguant par le fait, 
qu'on determine les valeurs seuils pour toutes les 
votes de haute frequence a I'aide de ('analyse de la 
repartition des valeurs des pixels de I'image premie- 
re. 

12. Le procede selon p. 1 se distlnguant par le fait, 
que le pixel c'est un vecteur, dont les composants 
caracter isent, par exempfe, les brillances des trois 
couleurs primitives (rouge, vert, bleu). 

13. Le procede selon p. 12 se dlstlnguant par le fait, 
qu'on calcu le lavaleurde la correlation pour chaque 
pixel par la multiplication scalaire de ce pixel a la 

. somme ponderee des pixels voisins. 

14. Le procede selon p. 13 et p. 5 se distlnguant par 
le fait, que pendant la calculation de fa somme pon- 
deree des pixels voisins on utilise les poids aniso- 
tropiques, la direction de fanisotropie correspon- 
dant a la direction de la transmission principale de 
la voie de frequence traitee. 



17. Le procede selon p. 1 se distlnguant par le fait, 
qu'avant de calculer des coefficients de pondera- 
tion on fait i'aplanissement des valeurs de la corre- 
lation pour plusieurs pixels voisins au moins dans 

5 une de (n-1) voies de haute frequence. 

18. Le procede selon p. 17 se distlnguant par le fait, 
qu'avant d'aplanir les valeurs de correlation on fait 
la transformation non-tineaire des correlations de 

10 telle maniere que les valeurs plus petites ou com- 
parables avec la valeur seuil restent inchangees, et 
les valeurs plus grandes diminuent. 

19. Le procede selon p. 1 se dlstlnguant par le fait, 
*5 que pendant ta formation du coefficient de ponde- 

ratton on fait I'aplanisseme de ce coefficient dans 
les pixels voisins au moins dans une de (n-1 ) voies 
de haute frequence. 

20 20. Le procede selon p. 1 se dlstlnguant par le fait, 
que I'image premiere represente une matrice des 
pixels de la dimension p, ou p £ 3. 

21. Le procede selon p. 1 se dlstlnguant par le fait, 
25 que pour la formation du coefficient de ponderation 
au moins dans une de (n-1 ) voies de haute frequen- 
ce on utilise les valeurs seuils differantes pour les 
parties de I'image differantes. 

30 22. Le procede selon p. 7 et p. 21 se distlnguant par 
le fait, qu'on determine les valeurs seuils pour les 
parties de ('image differantes et pour les voies de 
frequence differantes a I'aide de I'analyse de la re- 
partition des valeurs des pixels de I'image corres- 

35 pondants de la voie de frequence correspondante. 

23. Le procede selon p. 1 2 et p. 21 se distlnguant par 
le fait, qu'on determine les valeurs seuils pour les 
parties de I' image differantes et pour les voies de 
40 frequence differantes a ('aide de I'analyse de ta re- 
partition des valeurs ab so lues des vecteurs re pre - 
sentant les pixels de I'image correspodants de la 
voie de frequence correspondante. 

45 



15. Le procede selon p. 12 se dlstinguant par le fait, 
qu'on determine la valeur seuil pour chaque de (n- 
1 ) voies de haute frequence a I'aide de I'analyse de 
la repartition des valeu rs absotues des vecteurs, re- 50 
presentant les pixels de cette voie de frequence. 



16. Le procede selon p. 12 se distlnguant par te fait, 
qu'on determine les valeurs seuils pour toutes les 
voies de haute frequence a I'aide de I'analyse de la 55 
repartition des valeurs absolues des vecteurs re- 
presentant les pixels de I'image premiere. 
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